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Abstract 

Proposed low energy beta beam facilities would be capable of producing intense beams of neu- 
trinos (anti-neutrinos) with well defined spectra. We present analytic expressions and numerical 
results which accurately show how the total neutrino flux reaching the detector depends on the 
geometry of the source and the detector. Several authors have proposed measurements which re- 
quire using different flux shapes. We show that detectors of different sizes and shapes will receive 
neutrino fluxes with different spectral shapes, and that the spectral shape will also be different 
in different regions of the same detector. Our findings also show that for certain measurements 
systematic uncertainties and run time can be reduced. 
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I. INTRODUCTION 



The concept of a beta beam facility to produce neutrinos and anti-neutrinos has recently 
been proposed [l[ and feasibility studies have been performed [2'] and are underway . The 
source of the neutrinos will be radioactive ions which undergo beta decay. The ions will be 
boosted to relativistic speeds and contained in a storage ring which has two straight sections. 
Neutrinos emitted along a straight section will reach a detector placed along the axis of the 
straight section. The neutrino energy spectrum in the ion rest frame is known theoretically 
and verified by experiment. The spectrum from the boosted ions can thus be calculated and 
is well defined. The beta beam neutrino source is appealing because it will produce neutrinos 
of a single flavor [11]. Experiments to measure neutrino oscillation parameters, including the 
CP violating phase, where proposed in Refs 



A low energy beta beam was considered in [5| (for a recent review see Ref. [6|) to produce 
neutrinos with energy in the range of 10 to 100 MeV. This energy range is relevant for 
supernova neutrinos, and low energy neutrino physics in general. Typical boost factors that 
would generate neutrinos in this range are 5 < 7 < 15. Different sets of ring dimensions 
for a low energy beta beam, the resultant neutrino fluxes, and neutrino-nucleus inelastic 
scattering event rates at a close detector (~ 10 meters from the ring) were studied in detail 
in Ref. [7|. There it was found that the smallest possible ring, with as large as possible 
straight section for that ring, would result in a higher number of neutrinos reaching the 
target. This conclusion was independently confirmed in [8|. 

Several experiments to be done at a low energy beta beam have been proposed. These 
include measuring neutrino- nucleus inelastic scattering cross sections[3], testing the response 
of a lead based supernova neutrino detector]^, detecting a neutrino magnetic moment by 



measuring neutrino electron scattering [9|, measuring sin^iywith Vf. — e scattering [10|, test- 



ing conservation of the vector current by measuring v^. — 'p capture [ll(], reconstructing a 



detected supernova spectrum as a linear combination of beta-beam neutrino spectrum 12 
and measuring sin^^v^with neutrino- nucleus elastic scattering Tsl . A unique feature of the 
beta beam neutrino source is that it is possible to produce different spectra shapes with 
different energy ranges by boosting the parent ions to different 7 factors. This feature was 
first exploited in Ref. [sl and then later in References [lO, 12[ for the measurements they pro- 
pose. References 0, [l3| also suggest using this feature to measure cross sections at different 
energies. 

In this paper we derive, for a special case, an exact analytical expression that gives, for 
any point along the straight section, the percentage of the total emitted neutrinos that are 
contained in a lab angle Qj,. This expression is useful for considering detector sizes that will 
both receive a desired flux and maximize the number of neutrinos received. It is also useful 
for understanding the shape of the spectra at the detector. Using our expression, we give 
examples that show what percentage of neutrinos, from different points on the ring, reach 
the detector and what percentage of the neutrinos, from the same point on the ring, reach 
different parts of the detector. 

We consider for the first time the flux reaching detectors of various sizes. In particular, we 
calculate the neutrino flux reaching cylindrical detectors of the same mass but with different 
dimensions. We show that changing the radius and length of a cylindrical detector will 
alter the total flux and spectra shape reaching the detector. We also show that inner and 
outer regions of a cylindrical detector will receive different flux and spectra shape. Given this 
finding we propose that instead of running the beam at different energies in order to produce 
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a different shaped flux in the detector, the beam could be run at a flxed energy and events 
can be considered from different regions of the detector (that receive different flux shapes). 
This method to produce different shaped neutrino fluxes in a beta beam experiment has the 
advantage that no new systematic uncertainties are introduced from running the beam at 
different energies. This method could therefore allow for more accurate measurements and 
simplify fltting analysis. 

We calculate expected neutrino scattering events that would occur in the inner and outer 
regions of a cylindrical detector and compare these to events occurring from running the 
beam with ions boosted to different 7 factors. The event rates show that for certain mea- 
surements, more events can be obtained by using our method for producing different flux 
shapes in a detector. Thus our proposed method for producing different shaped neutrino 
fluxes not only can reduce systematic error, it can also reduce the beam time needed. 

The average energy of neutrinos emitted from a beam of relativistic particles will have 
some dependence on the angle of emission relative to the beam axis. This is a consequence 
of special relativity, and we take advantage of it in our considerations of the neutrino fluxes 
reaching different regions of an on-axis detector. This property can also be exploited for use 



in another context — putting a detector off a beam axis, e.g. for pion decay in flight |14j . |15 
and recently for a beta beam 161 . 

The paper is organized as follows. In section II we derive the expression for the percentage 
of neutrinos emitted from a point on the straight section that are contained in a given lab 
angle. We also present flgures to demonstrate the use of this expression. In section III 
we discuss different detector shapes and present plots of the calculated flux reaching these 
detectors. In section IV we examine the flux reaching different regions of a detector and we 
calculate neutrino scattering event rates in these regions. We show how certain experiments 
can be made more efficient by exploiting our findings regarding the fiux reaching seperate 
regions of the detector. We give conclusions in section V. Our formalism for calculating the 
neutrino fiux in the lab frame is presented in the appendix. Throughout the paper we use 
units in which h = c = 1. 



II. NUMBER OF NEUTRINOS EMITTED IN A LAB ANGLE 

A currently used estimate for the amount of emitted fiux in the lab frame is that most of 
the fiux is contained in an angle of I/7, where 7 is the boost factor of the ions. For example, 
this estimate was used in Ref. 0]. (Note that Ref. [7] also provided exact calculations of 
numbers of events.) The estimate may be used to choose detector size parameters suitable 
for capturing most of the fiux coming from the straight section. However, one may desire a 
more accurate estimate. In particular, we will show that maximizing the number of neutrinos 
that reach the detector will not necessarily maximize the fiux. 

We aim to quantify how many neutrinos are emitted in a lab angle. More specifically, 
given any point along the straight section of the storage ring, what fraction of the total 
neutrinos emitted from that point at any instant in time are contained in a given lab angle? 
In the appendix, we present our formalism which expresses the probability that a neutrino 
will be emitted with a certain energy and angle in the lab frame. To answer our present 
question we refer to Eq. (1A6P and the coordinate transformation of Eq.s (lAlip - (1A13P and 
consider a special case. Instead of using the rest frame beta decay spectrum of Eq. (lA6p . 
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we will use a monoenergetic rest frame spectrum defined by a delta function 



fR{En)=C6{En-Eo),. (1) 

Here, C has dimensions of per second and the delta function has implicit dimensions of per 
energy. Using Eq. (lASP and transforming this flux in the manner of Eq. (1A17I) we have 



1 /■^r'' C S[EL-f{l - vcos Ol) - Eq] dEL dcos 9l 
-Jo 2C 7(1 - vcosOl) 



(2) 



where v is the relativistic speed of the nuclei and 7=1/ ^/l — v'^ is the boost factor. We know 
from the discussion in the appendix Sec. IA2I that if we restrict the bounds of integration 
in Eq. ([2]), then Eq. ([2]) will give normalized probability that a neutrino will be emitted 
with lab energy G {El,El) and lab direction cosine G (cos cos 6*1). For the case of a 
monoenergetic spectrum there is only one energy represented in the rest frame. However, 
since the rest frame angle is arbitrary, that one energy can be transformed to a range of 
energies in the lab frame. Equation ([2]) is useful since the integrations can be performed 
analytically using the identity 



df 




dy 





5{y-yi), (3) 



where yi are simple zeros of f{y)- 

The upper and lower bounds on the cos 6' integration in Eq. ([2]) correspond to angles 
6 = and 6 = n, respectively. We may change the lower bound to cos ^"P'^'^ corresponding to 
some opening angle 6'°^'^'^. In this case, performing the integrations in Eq. ([2]) would give the 
probability that a neutrino is emitted in the lab frame between 6'°p™ and ^ = 0. Denoting 
this probability P°p'^^^ we have 

popen ^ I -COS 6^^- 

272(1 COS e°P°°)' ^ ^ 

Notice that this result depends only on the speed v to which the ions are boosted. This 
makes sense given Eq. (1A14|) which shows that the lab frame angle only depends on v and 
cos Or, the angle the neutrino was emitted from in the rest frame. If the neutrino is emitted 
in a more forward direction in the rest frame {cos Opt closer to 1), Eq. (lAlip implies the 
neutrino will have a higher energy in the lab frame. Thus, the amount of neutrinos that get 
forward boosted within a given lab angle depends on the speed of the nuclei, and neutrinos 
with smaller lab angles will have a higher lab energy than do neutrinos with larger lab angles. 
These considerations apply to our current case of a monoenergetic rest frame spectrum, but 
they generalize to the case of the distribution spectrum given by Eq. flA6p . In the next 
section we will see the consequences, for the distribution spectrum, of the fact that more 
forward boosted neutrinos have higher energies than do neutrinos with large lab angles. 

In the appendix Sec. lA 41 we explain that uj/L (where uj is the number of nuclei/sec 
injected in the storage ring and L is the total length of the ring) is the number of neutrinos 
emitted per second from a line element dx on the ring. Multiplying P^p'^^ hj u/L gives, for 
any point, the number of neutrinos per second emitted inside angle ^"P'^'^, where the angle 
is measured from the point of emmission. We now define ^°p™, the fraction of neutrinos 
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emitted from some point on the straight section that are contained in angle ^"p*^"^. Since the 
total number of neutrinos per second emitted from any point on the straight section is u/L 
we have that 

o^popen 

So we see that ^"p*^" = p°pen qq^^i now calculate, for the case of a monoenergetic 

rest frame spectrum, exactly what fraction of the neutrinos emitted from any point on the 
straight section are contained in a lab angle 6l ~ I/7 for a beam run at 7 = 7. Using ^"p*^^ 
we find that only 50% of the neutrinos emitted at each point are contained in this angle for 
this case. 

To illustrate the use of ^"p*^'^, we consider two examples of detector configurations. The 
dimensions of the detectors in these examples are not intended to suggest an actual experi- 
mental setup, but rather to demonstrate the use of ^"p^'^. In the first example we show the 
fraction of neutrinos emitted from two different points on the straight section that reach 
the same region of the detector. In the second example we show the fraction of neutrinos 
emitted from the same point on the straight section that reach two different regions of the 
detector. These examples are shown in Fig. ([2]) and Fig. ([3]), which are to scale. In each 
figure, the thick black line is the straight section of the storage ring, the rectangle is a side 
view of the cylindrical detector, and the dashed cone shows the opening angle that contains 
90% of the neutrinos emitted from that point. 

In Fig. ([2]), we show a detector with radius 6m and thickness 0.5m. The detector is 10m 
away from the straight section, which itself is 50m. The opening angle from point B to the 
detector is 16.7° which, using ^^p^^"^^ contains 80% of the neutrinos emitted from that point. 
Thus, only 80% of the neutrinos emitted from that point reach the detector. The opening 
angle from point A to the detector is 6.24°, which contains 37% of the neutrinos emitted 
from that point. 

In Fig. ([3]), we show a detector with radius 6m and thickness 14m that is located 10m 
away from the straight section. The opening angle from point C on the straight section to 
the detector slice at D is 14°. Using ^°p'='^^ we calculate that 75% of the neutrinos emitted 
from point C reach this slice of the detector. The angle from point C to the detector slice 
at E is 10° and 60% of the neutrinos emitted from point C reach this slice. This figure also 
shows that the neutrinos emitted between opening angle of 14° and 10° only have a small 
portion of the detector to interact with. 



III. MODIFYING DETECTOR DIMENSIONS 
A. Number of Neutrinos vs. Flux 

In the previous section we saw how to quantify the number of neutrinos reaching the 
detector. We now study how changing detector dimensions will effect the neutrino flux in 
the detector. It is apparent from Fig. ([2]) that a larger radius will allow the detector to 
catch more neutrinos. However, the flux is the number of neutrinos (per second per energy) 
reaching a slice of the detector divided by the area of that slice. Increasing the radius has 
two effects. Allowing more neutrinos to reach the detector works towards increasing the 
flux. Dividing by a larger area will reduce the average flux. 

It is apparent from Fig. ([3]) that if the thickness of the detector is increased less neutrinos 
will reach the back end of the detector than will reach the front end of the detector. However, 
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the neutrinos that do reach the back of the detector will have had opportunity to interact 
with more detector material than do neutrinos emitted with larger lab angles. The neutrinos 
that are exposed to more detector material will have a greater chance of having an interaction 
within the detector. 

In the case of having a fixed thickness, increasing the radius of the detector will increase 
its total mass. Likewise, for a fixed radius, increasing the thickness will also increase the total 
mass of the detector. Larger mass detectors of course have more targets and the potential 
for more events. 

Motivated by Figures ([21)- © and these qualitative considerations, we use Eq. flA27p to 
calculate the average flux in detectors of different radii and thickness. We would like to 
study whether a long detector with small radius (see detector 1 in Fig. (jl])) would be more 
desirable than a short detector with a large radius (detector 2 in Fig. (j4])). When comparing 
detectors with different radii and thickness it will be useful to have some common property 
of the detectors. Therefore, we will keep the volume fixed and vary the radius and thickness 
while maintaining the volume. For our example, we chose a detector with volume 71.2 m^ 
and we use five different pairs of radius r and thickness h, both in meters. The pairs we 
have chosen, denoted as {r,h), are (2.13, 5), (2.75, 3), (3.89, 1.5), (4.77, 1) and (6.74, 0.5). 
We calculate the fluxes for the case of a beta beam facility with storage ring of total length 
L = 450 m, straight section length S = 150 m, detector located d = 10 m away from the 
straight section and with a; = 5 x 10^^ nuclei/sec injected into the ring 0]. We do our 
calculations for a beam composed of ^^Ne nuclei. ^^Ne has a half life of 1.67 seconds and 
I3~^ decays emitting a neutrino with a maximum possible energy of 3.42 MeV. We use this 
configuration of the beam for all our calculations throughout the paper. The fluxes reaching 
the detectors mentioned above are plotted in Fig. (jSj) for the cases of a beam run at boost 
of 7 = 7 and 7 = 14. For comparison, in Fig. ([6]) we show different flux plots calculated for 
a beam run at boost factors of 7 = 7, 8, 9, 12, 14 with a 71.2 m^ detector with r=2.13m and 
h=5m. 



B. Flux Curves 

The plots in Fig. ([6]) demonstrate the attractive feature of the beta-beam neutrino 
experiment, that it is capable of producing neutrino fluxes of different shapes by boosting 
the nuclei to different 7 factors. This feature was utilized for the measurements proposed 
in References js], [l^, 12]. References 0, 13] also mentioned taking advantage of this feature. 



The plots in Fig. ([5]) illustrate our conclusion that another way to produce different shaped 
neutrino fluxes is by using different shaped detectors with a beam run at only one boost 
factor. 

From Fig. (jSj) we see that for the detector dimensions we have chosen in this example, the 
flux decreases as the radius of the detector increases. This is consistent with our qualitative 
observations above. When considering detector sizes, beyond technical and engineering 
restrictions, one should strive to find a balance between maximizing the number of neutrinos 
reaching the detector and maximizing the flux within the detector. 

We also notice from Fig. ([5]) that as the radius decreases the average energy for the flux 
distributions decreases. This is consistent with what we found in the previous section for the 
case of a monoenergetic rest frame spectrum - that neutrinos emitted with smaller angles in 
the lab frame, i.e. more forward boosted, have higher energies. Obviously, the distribution 
rest frame spectrum is more complicated and does not behave exactly as a mono-energetic 
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spectrum does when transformed to the lab frame. However, Fig. ([5]) clearly shows that 
there are more higher energy neutrinos in the forward boosted direction. 

We see from Fig. that if we produce different flux distributions in the detector by 
changing the size of the detector, we will always have neutrinos with energy in the same range 
[0, EJ^^^]. Comparing to Fig. where different flux distributions are produced by boosting 
the nuclei to different 7 factors, we see that the available energy range of the neutrinos 
decreases as 7 decreases. Thus, changing the spectrum by changing the detector could be 
useful if we want to compare different measurements involving higher energy neutrinos. 

Another potential advantage of changing the flux shape by changing detector dimensions 
is that no new beam systematic uncertainties are introduced. If the beam is run at different 
boost factors there may be different losses during acceleration, or different losses due to 
nuclei decaying before reaching the ring. However, if uncertainties stemming from using 
different sized detectors are lower than beam uncertainties, it would be favorable to produce 
different shaped spectra by running the beam at a fixed energy and using different sized 
detectors. 

IV. FLUX SHAPES AND EVENT RATES IN DIFFERENT REGIONS OF DE- 
TECTOR 

A. Flux Shape in Separate Regions of Detector 

We now present another new way for obtaining different flux shapes in a beta beam 
experiment. Building on the results of the previous two sections, we investigate the shape 
of the flux in different regions of one detector. We divide a cylindrical detector into inner 
cylinder and outer hollow cylinder regions, as in Fig. ([7]), and find that different flux shapes 
reach these seperate regions. This may not come as a surprise given the results shown in the 
previous two sections. However, to illustrate this result we present plots of calculated fluxes 
for the inner and outer regions of a cylinder detector. Continuing with our current example, 
we use a 71.2 m^ detector with r=2.13m and h=5m. We divide this detector into seperate 
inner cylinder and outer hollow cylinder regions. We consider two cases for separation: 
one with the boundary at radius r=lm and another with the boundary at radius r=1.5m. 
For the case with boundary at r=lm, the inner cylinder has volume 15.7 m^ and the outer 
hollow cylinder region has volume 55.6 m'^. For the case with the boundary at r=1.5m, each 
region has a volume of about 35.5 m'^. In each case, we use the same parameters for the 
experimental setup that we have been using: storage ring of total length L = 450 m, straight 
section length S = 150 m, detector located = 10 m away from the straight section and 
a; = 5 X 10^^ nuclei/sec injected into the ring. We show plots of the calculated fluxes for 
these two cases, for nuclei boosted to 7 = 14, in Fig. ([8]). 

We can see from Fig. ([8]) that the neutrino flux in the inner and outer regions of a detector 
do have different shapes. In each figure we see that the flux in the inner region has a higher 
peak than the flux over the whole detector. Thus we see that there is a greater concentration 
of the neutrinos in the inner region of the detector. We see also that the average energy of 
neutrinos reaching the inner region is greater than the average energy of neutrinos reaching 
the outer region. We therefore conclude that another way to achieve different shaped fluxes 
in a beta beam experiment is to separately count events that occurred in an inner region 
and outer region of a detector. Location for events can be identified in Cherenkov and 
scintillator detectors. Instead of running the beam for one year at some boost factor and 
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then another year at a different boost factor, or running the beam for one year with one 
detector and another year with a different sized detector, the beam could be run for two 
years with one detector and events could be taken from two seperate regions. This method 
has the advantage that no new systematic uncertainties are introduced from using different 
detectors or different beams. 



B. Events in Seperate Regions of Detector 

We illustrate this method for obtaining different flux shapes by calculating event rates 
for neutrino-nucleus elastic scattering. This process was studied in Ref. at a beta beam 
facility and also in Ref. [13] for neutrinos produced at a stopped-pion source. For our exam- 
ple we consider a liquid noble nuclear recoil detector based on CLEAN [18]. This detector 
concept has been proposed with a spherical geometry, but cylindrical geometries would also 
be possible [3]. We calculate events for a 100 tonne fiducial volume in a cylindrical detector 
containing liquid argon. Liquid argon has a density of 1.4g/cm^ so 100 tonnes will occupy 
a volume of 71 m^. For our example we will choose a fiducial region with radius of 2.13m 
and thickness of 5m. An actual detector would have a larger total mass and larger total 
dimensions and the region surrounding the fiducial volume would be used to reject some 
background events from external sources. Since the purpose of our example here is to illus- 
trate how different flux shapes could be achieved by separating events from different regions 
of a detector, and not intended for proposing a specific measurement, we do not address 
the issue of backgrounds that would exist in the detector we are considering. Our finding 
is general and applies not only to liquid noble detectors, but any detector that has position 
resolution for neutrino interaction events. 

This method applies to detectors of other shapes as well. For example, recalling from 
Section lTlIBI that forward boosted neutrinos have higher energy, a cone shaped detector (such 
as detector 3 in Fig. (jlj)) would be useful for selecting higher energy neutrinos. Engineering 
and technical limitations may prevent a realistic detector of this shape but such a detector 
could be simulated in a cylindrical detector. By selecting events that occur within a cone 
shaped region of a detector we could pick out events which involved mainly higher energy 
neutrinos. 

Returning to the cylindrical detector, we calculate the number of neutrino-nucleus elastic 
scattering events per year per nuclear recoil energy by folding the average flux (Equation 
IA27|) with the differential cross section for neutrino-nucleus elastic scattering: 



1 71 r ^max i 



Here, E is the neutrino energy in the lab frame, T is the nucleus recoil energy, min(T) = 
(T + VWT2TM)/2 is the minimum neutrino energy required to give a nucleus of mass M 
recoil energy T, and Nt is the number of targets in the detector. We use the cross section 
given in Ref. [19|]. We neglect radiative corrections and use the analytic form factor given 
m Ref. 0. Calculated event rate curves are given in Fig. ([9]) and Fig. (fTOj) . In Fig. ([9]) we 



show events in the inner cylinder and outer hollow cylinder regions for a fiducial region of 
100 tonnes. The boundary of the regions for the plots in the left panel of Fig. (Q is at r=lm; 
the inner and outer regions have masses of 22 and 78 tonnes, respectively. The boundary 
of the regions in the right panel of Fig. (Q is at r=1.5m; the inner and outer regions each 
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have mass of 50 tonnes. For comparison, we show in Fig. (fTOj) the events occurring in the 
entire fiducial volume for the case of a beam run at two different boost factors, 7 = 7, 9. 

We see from Fig. (Q that event rate curves will be different for events occurring in 
different regions of a detector. Notice that for the case of the boundary at r=lm, the 
left panel of Fig. there are more events in the outer region, while for the case of the 
boundary at r=1.5m, the right panel of Fig. there are more events in the inner region. 
The average flux in the inner region is highest for the case of the boundary at r=lm, but this 
region has the smallest mass and hence less targets. This is consistent with our qualitative 
observations in the previous section. 

Recall the importance of different flux shapes for the measurements proposed in Refer- 
ences 0, S, [13, 12, 13|. To illustrate our method, suppose a measurement required a beam 



to be run for a year at 7 = 13 and then another year at 7 = 14. The measurement could 
likely also be accomplished by running the beam at 7 = 14 using one detector and taking 
events from two separate regions for a duration of two years. We now discuss how using 
our method to produce different flux curves in a detector could also reduce total run time 
for some experiments. Suppose a measurement required running the beam for a year at 
7 = 7 and 7 = 9. Instead, the beam could be run for one year at 7 = 14 and events could 
be taken from different regions of the detector. This situation is demonstrated by flgures 
(Q and (ITO!) . Each curve in Fig. (ITOl) would require a year of beam time - for a total of 
two years for the measurement. However, only one year is required to produce both dashed 
curves in, for example, the right panel of Fig. (Q. The dashed curves in the right panel 
of Fig. not only show more events, they show events for higher recoil energies as well. 
Thus we see with our method, not only can systematic uncertainties be reduced, for certain 
measurements more events can be obtained in less time. 



V. CONCLUSIONS 

We have derived an expression which offers a way to quantify the fraction of neutrinos 
emitted from a point on the straight section that are contained in a given angle in the lab 
frame. We have calculated flux shapes and events rates for different detector conflgurations. 
Different flux shapes can be obtained by running the beam at a flxed energy and using 
detectors with different dimensions, or by using different regions of a single detector. The 
latter method has some advantages, including providing more events for certain measure- 
ments and reduced uncertainties. Although detailed studies of uncertainties remain to be 
done, we anticipate that the latter method will have less systematic uncertainties than the 
methods of running the beam at different energies or using different shaped detectors. 
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APPENDIX A: CALCULATING NEUTRINO FLUX IN LAB FRAME 



1. Rest Frame Spectrum 



We start with the beta decay rate in the nucleus rest frame [21 



^=1 -^TT7:F{±Z,E,)E,p,EldE,. (Al) 

Here rrie is the electron mass, {ft) is the ft-value containing the nuclear matrix element, E^ 
and Pe are the electron energy and momentum respectively, E^, is the neutrino energy (we set 
the neutrino mass mj, = 0), and F{±Z,E(,) is the coulomb correction factor. The coulomb 
correction factor accounts for the the coulomb interaction between the emitted electron (or 
positron) with the charge of the nucleus. 

The energy of the emitted electron depends on the energy of the emitted neutrino, 

Ee = Qnuc — E^, (A2) 

where Qnuc is the nuclear mass difference (not to be confused with the atomic mass difference) 
for the transition. The electron momentum is determined from the electron energy by the 
usual relativistic expression 



p, = ^El-ml. (A3) 

The neutrino has maximum energy when the electron has zero momentum. From Eq. flA2p 
we see this is 

ET"" = Qnuc - me. (A4) 

The half life for the decay is 

ln2 



ti/2 = — = rln2, (A5) 



where r = 1/A is the lifetime. 
The integrand of Eq. ( lAll) is 



ln2 

fR{Eu) = -^rj-F{±Z, E,)E,p,El (A6) 

We will refer to the function ffi{Ey) as the neutrino energy spectrum in the rest frame. Note 
that 

— — ~ ~ = -TfRiEu)dEi, (A7) 

fR{E,)dE, 



is the normalized probability that the neutrino will be emitted with energy in the range 
{Eu, Ey + dEy) when a nucleus undergoes /3-decay. The expressions on each side of Eq. ( ]A7p 
are dimensionless. 

There is no preferred direction for the neutrino to be emitted in the nucleus rest frame; 
the decay is isotropic in direction. For the coordinate system defined in Fig. ([1]), ^ G [0, tt] 
and if G [0,27r], there is an equal probability of the neutrino being emitted in the range 
(0, 9 + dO) for any angle 9. We can define a new function for the neutrino spectrum in the 
rest frame which is also a function of the direction cosine coordinate cos^^ G [—1, 1]: 

fR{E,,cos9) = fR{E,). (A8) 
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This function is constant in cos^. Note that 

/ / fR{E^,cose)dE^dcose = 2X. (A9) 
J-i Jo 

Now we have 

l^^^^^^dK,dcose (AlO) 
2 A 

for the normahzed dimensionless probabihty that the neutrino is emitted with direction 
cosine G (cos 6*, cos 6* + dcos9) and with energy e {E^, E^, + dE^). 



2. Transformation to Lab Frame 

Now consider the rest frame moving with velocity v with respect to the lab frame. Re- 
ferring again to figure Fig. ([1]), we take the nucleus rest frame to be moving in the direction 
of the arrow on the axis shown. The neutrino emitted in the rest frame has energy E^ and 
the longitudinal and transverse components of momentum are and p^, respectively. (We 
have dropped the subscript v for clarity.) The lab frame energy and momentum components 
are: 

El = -fEnil + VCOS 9r) (All) 
pi = lEniy + cos^/j) = ElcosOl (A12) 
pi = E lysine R = ELsm6L. (A13) 

Recall that v = v^7^ — I/7. Equations IA12I and lAllI can be combined to obtain the expres- 
sion 

cos6'l = {cos 9r + v)/{l + vcos 9r). (A14) 

To derive an expression for the neutrino spectrum in the lab frame, we start with Eq. 
(]A10|) . The transformation equations (lAlip - (]A13I) can be inverted to give the rest frame 



coordinates in terms of the lab frame coordinates: 



Er = Er{El, cos Ol) (A15) 
cos9r = cos9r{cos9l)- (A16) 



Then we have 

fniER, cos 9 R 



2 ^ dER rfcos 9r — > 
fR[ER{EL, cos 9l), cos 9r{cos 9l)] O^Er, cos 9r] 



dELdcos9L, (A17) 



2 A d{EL,cos9L) 
where the Jacobian is d{ER, cos 9R)/d{EL, cos 9 l) = 1/7(1 — vcos9l)- Finally we have 

f fj? a \AT7 A a fnlERiEL, cos 9l)] dEi dcos 9l , . 

fL{EL,cos9L)dELdcos9L = — J-. T-^. (A18) 

2 A 7(1 — vcos 9l) 
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For an example to demonstrate this coordinate transformation we check the normahzation: 



= t r^^^ ^-^^dEndcosOn (A19) 

'■^ r^r"- fR[ER{EL,COS0L)] dEi^dcOS0L . , 



1 JO 



2 A 7(1 — vcos^l) 

I'll''' fL{EL:Cos9L)dELdcos9L. (A21) 



-1 ^0 

For a /?-dccay nucleus boosted to speed v we interpret //.(-El, cos^L)(i£'L (icos^L as the 
normalized probability that the neutrino will be emitted with lab energy G {Ei, Ei + dE^) 
and lab direction cosine G (cos 6*^, cos 6'^ + dcos9L)- We illustrate this with a thought 
experiment. Consider a large spherical neutrino detector. We send a nucleus boosted to 
speed V through a central axis of the detector. The detector is large enough that the nucleus 
will decay before it reaches the end of the axis. We are able to identify the position of 
the nucleus when it decays and we detect the neutrino in the detector. Thus we know the 
lab direction cosine that the neutrino was emitted with. We also measure the neutrino's 
energy when we detect it. Now we repeat the experiment with another nucleus over and 
over. The distribution of the "measured" neutrino energy and direction cosines will be given 
by f l{E L-, cos 9 l). 



3. Number of /3-decays per second in the Ring 

Consider the storage ring. Initially, at time t — there are zero nuclei in the ring. We 
start injecting uj ions per second into the ring at i = 0. We would like to know the total 
number of decays per second that occur in the ring after an equilibrium situation has been 
reached. In the lab frame, the decay rate is Al = l/r^, where tl is obtained from the 
lifetime in the rest frame by tl = 'JTr. Let N{t) be the number of nuclei in the ring at time 
t, and dN be the change in the number of nuclei during time interval dt. Then we have 

dN = -NXLdt + ujdt. (A22) 

Using the initial conditions, the solution to this differential equation is 

N{t) = ^ - ^e-^^K (A23) 
Al Al 

Equilibrium is reached in the limit t — > 00. The number of nuclei in the ring, N{t), after 
equilibrium is reached is 

lim N(t) ITruj. (A24) 

t^co Al 

The number of decays in the ring, XLN{t), after equilibrium is reached is 

lim XlNU) = u. (A25) 

t— >oo 

After a long time of injecting ions into the ring, the number of /3-decays occurring in the 
entire ring is equal to the number of ions per second injected into the ring. 
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4. Integrating Over the Ring and Detector 



A cylindrical detector, with radius r and length is placed at a distance d away from the 
ring along the axis of one of the straight sections. The flux at each point in the detector is 
not the same since neutrinos are incident from different angles. To compute the number of 
events in a detector we need to fold the flux with the cross section for the neutrino scattering 
interaction. We could calculate the flux at each point in the detector, multiply the flux at 
each point by the cross section, integrate over position in the detector, and then integrate 
over neutrino energy. However, we are not interested in any angular dependence of events. 
Rather than calculating a flux at each point in the detector, we will use the average of the 
flux at all points. To get the total events we take the average flux at a point, multiply by 
the number of targets in the detector (i.e. the number of points), and fold this with the 
cross section. For the purposes of calculating total events this is equivalent to the preceding 
way to do the calculation. However, one should keep in mind that if angular dependence of 
the events is desired, perhaps for considering detection mechanisms, the former method is 
necessary. 

We first consider the case of an infinitesimally thin detector. To calculate the average flux 
in this detector we add the total number of neutrinos per second passing through the slice 
and then divide by the area of the slice. Let N^{El) (IEl be the total number of neutrinos 
per second, with lab energy G {El, El + cIEl), passing through detector. 

The /3-decays occur throughout the entire storage ring. For a ring of total length L, 
with straight section lengths S, the number of decays per second per length element dx 
occurring at any point along the ring is uj/L. To find Nu{El) dEL we sum the contribution 
of neutrinos from points along the straight section, (u/L) /l^El, cos 9L)dELdcos9L, that 
reach the detector: 



N,{EL)dEL = u / fL{EL,coseL)—dcoseLdEL. (A26) 

J-=^±±^ Jo L 



For a detector with thickness h, Eq. (]A26I) (with x + d x + d + z) is the number of 
neutrinos per second passing through the slice of the detector located between {z,z + dz). 
Integrating over dz/h gives the average of these numbers over all slices. Finally, the average 
of the flux over all points in the detector is obtained by dividing this average over slices by 
the cross sectional area of the detector: 

^El) = —J' [' ['fLiEL,coseL)^^dcoseL. (A27) 

7rr^ ' , ^+'+'1 Jo Jo h L 

Notice that Eq. flA27l) has dimensions of number per energy per time per area. 
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FIG. 1: Coordinate system defined in the rest frame of the nucleus. 
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FIG. 2: (Color Online) Example configuration to show the use of to calculate the fraction 
of neutrinos emitted from two arbitrary points along the straight section, that reach the detector. 
This example is for nuclei boosted to 7 = 7. The figure is to scale. The thick black line represents 
the straight section, the rectangle is a side view of the cylindrical detector. The dashed cones show 
the lab angle that contains 90% of the emitted neutrinos. The dashed-dot lines show the opening 
angle to the detector face. 




FIG. 3: (Color Online) Example configuration to show the use of ^"P'^'^ to calculate the fraction 
of neutrinos emitted from an arbitrary point along the straight section, that reach two arbitrary 
slices of the detector. This example is for nuclei boosted to 7 = 7. The figure is to scale. The thick 
black line represents the straight section, the rectangle is a side view of the cylindrical detector. 
The dashed cones show the lab angle that contains 90% of the emitted neutrinos. The dashed-dot 
lines show the opening angle to the two detector slices. 
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FIG. 4: Examples of shapes for a detector. 



Neutrino Flux from y = 7 Beam in 71 .2m^ Detector Neutrino Flux from y= 14 Beam in 71 .2m^ Detector 
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FIG. 5: (Color Online) Neutrino fluxes in detector with volume 71.2 m^ for five pairs of cylinder 
radius r and thickness h, both in meters. Pairs are denoted as (r, h) in the key. Fluxes were 
calculated for storage ring of total length L = 450 m, straight section length S = 150 m, detector 
located at d = 10 m away from the straight section and a; = 5 x 10^^ nuclei/sec injected in the 
ring. In the left panel, fluxes are calculated for a beam run at a boost of 7 = 7; in the right panel, 
fluxes are calculated for a beam run at a boost of 7 = 14. 
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FIG. 6: (Color Online) Neutrino fluxes, for beam run at boost factors of 7 = 7, 8, 9, 12, 14, in 
detector a 71.2 m^ detector with r=2.13m and h=5m. Experimental setup assumed is the same as 
that for the fluxes calculated in Fig. ([5]). 
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FIG. 7: Cylinder detector separated into an inner cylinder and an outer hollow cylinder region. 
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Vg Flux in (r=2.13m, h=5m) Detector, boundary r=1.5m 
10000 



E 
o 

> 



o 



E 



10 20 30 40 50 60 70 80 90 100 
Neutrino Energy (MeV) 




10 20 30 40 50 60 70 80 90 100 
Neutrino Energy (MeV) 



FIG. 8: (Color Online) Neutrino fluxes in detector with volume 71.2 and r=2.13m and h=5m. 
Fluxes were calculated for separate inner and outer regions of the detector. Pairs in the key denote 
the radii of the cylinder regions - the first number is the location of the outer radius and the second 
number is the location of the inner radius. Fluxes were calculated for a storage ring of total length 
L = 450 m, straight section length S = 150 m, detector located at d = 10 m away from the straight 
section, w = 5 x 10^^ nuclei/sec injected in the ring and a beam run at a boost of 7 = 14. The solid 
curve shows the flux in the whole region, the long dash curve shows the flux in the inner region 
and the short dash curve shows the flux in the outer region. For the plots on the left panel, the 
boundary between the inner and outer regions is located at r=lm; for the plots in the right panel 
the boundary is at r=1.5m. 
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FIG. 9: (Color Online) Number of neutrino-nucleus elastic scattering events per year per nuclear 
recoil energy in 100 tonne fiducial volume of liquid Ar. Pairs in the key denote the radii of 
the fiducial volume regions - the first number is the location of the outer radius and the second 
number is the location of the inner radius. Beam configuration is the same as for other calculations 
throughout the paper, with beam run at 7 = 14. Solid line shows calculated events in entire fiducial 
volume, long dashed line shows events in inner region and short dashed line shows events in outer 
region. Left panel is the case with boundary between inner and outer region at r=lm, inner region 
containing 22 tonnes Ar and outer region containing 78 tonnes. Right panel is for the case with 
boundary at r=1.5m, each region containing 50 tonnes of Ar. 
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FIG. 10: (Color Online) Number of neutrino-nucleus elastic scattering events per year per nuclear 
recoil energy in 100 tonne fiducial volume of liquid Ar. Experimental configuration is the same as 
for other calculations throughout the paper. Solid line shows calculated events for beam run at 
7 = 7, long dashed line shows events for beam run at 7 = 9. 
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